How Flow Changes Polymer Depletion in a Slit 
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Abstract. A theoretical model is developed for predicting dynamic polymer depletion effects under the 
influence of fluid flow. The results are established by combining the two-fluid model and the self-consistent 
field theory. We consider a uniform fluid flow across a slit containing a solution with polymer chains. The 
two parallel and infinitely long walls are permeable to solvent only and the polymers do not adsorb to 
these walls. For a weak flow and a narrow slit in ©-solvent conditions, an analytic expression is derived to 
describe the steady state polymer concentration profiles. In both 0- and good-solvents, we compute the 
time evolution of the concentration profiles for various flow rates characterized by the Peclet number. The 
model reveals the interplay of depletion, solvent condition, slit width, and relative strength of the fluid 
flow. 



1 Introduction 

In a polymer solution near an interface the polymer seg- 
ments are either attracted or repelled by that interface 
[1]. In the latter case there exists a depletion zone near 
the interface. In this zone the polymer segment concen- 
tration is smaller than the bulk value because the possi- 
ble number of configurations of the polymer chains is re- 
duced. The non-adsorbing wall forbids a certain amount of 
paths of the polymer chain. For ideal polymer chains near 
a hard wall the depletion thickness is close to the poly- 
mer's radius of gyration [2]. This result holds generally for 
a dilute polymer solution [3,4], whereas in a semi-dilute 
polymer solution the depletion thickness is determined by 
the correlation length [5] . These results for the dilute and 
semi-dilute concentration regimes were later combined [6, 
7]. The investigations of polymer depletion at an inter- 
face were primarily focused on the equilibrium case in the 
past. The equilibrium depletion thickness suffices to de- 
scribe the attraction between colloidal particles when the 
depletion layers overlap [8,9], which can be measured for 
instance by optical tweezers [10]. The depletion force may 
yield phase transitions [11-13] for which the binodals can 
be predicted for well-defined colloid-polymer mixtures [14, 
4]. It is of fundamental and practical interest to under- 
stand the change of the depletion layer under a fluid flow 
effect (Fig.l). Simple shear flow of a polymer solution next 
to a single wall leads to slip even if the depletion layer is 
assumed to be unaffected by shear [15]. It has been shown, 
however, that the flow does affect the depletion thickness 
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Fig. 1. Convective depletion effect for polymer chains near a 
non-adsorbing solid surface (a) under no flow and (b) under 
flow, n and t denote unit normal and tangential vectors to the 
surface. 



at a wall [24] . At the particle level, the models developed 
for the single particle motion and pairwise particle interac- 
tions [25-27] neglect the slight distortion of the depletion 
zone, which is applicable for long-time Brownian diffusion 
with weak convective effect. A convective depletion model 
was first established by Odijk [28] for describing a thin de- 
pletion boundary layer in front of a fast moving sphere. A 
complete picture of the fully coupled convective depletion 
effect under various solution conditions is not yet avail- 
able. 

The flow of a polymer solution near an interface is an 
important practical phenomenon [16, 17]. To give a few ex- 
amples, in oil recovery, aqueous polymer solutions are used 
to sweep oil out of natural porous media in the earths' sub- 
surface [18]. In that situation the polymer solution flows 
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through narrow pores where a large effective viscosity, as 
enhanced by the presence of the polymers, is required to 
flush oil out. Measurements of apparent polymer solution 
viscosity obtained from flow through narrow pores indi- 
cate that the viscosity is reduced compared to the bulk 
viscosity, see e.g. Chauveteau [19] and Omari et al. [20]. 
Capillary rheometry is another example [21]. It uses the 
fact that the flow time is proportional to the liquid viscos- 
ity. However, for polymer solutions the relative viscosity 
from capillary rheometry may differ from other methods 
[21-23]. These are all related to dynamic depiction effects 
in confined geometries. 

In this paper we propose a theoretical framework to 
investigate dynamic depletion effects. We study the influ- 
ence of flow on the polymer segment density profile in a 
narrow or wide slit. This is an interesting problem since for 
a narrow slit an analytical expression is available only for 
the equilibrium segment density profile [29,1,6]. It is un- 
clear yet to what extent a flow field modifies the depletion 
layers. To resolve the interplay of convective and diffusive 
effects, we start from a general formulation of the two- 
fluid model proposed by Doi and Onuki [30] for polymer 
transport. Then we apply the ground state approximation 
(GSA) to the self-consistent field theory (SCFT) and the 
dynamics version of SCFT to express the chemical poten- 
tial of polymer segments under various flow conditions. 
This framework allows us to clarify how the segment den- 
sity is altered upon increasing flow. It will tell for which 
Pcclct numbers fluid flow significantly affects polymer de- 
pletion as compared to the equilibrium case. A detailed 
analysis for 0- and good-solvent cases is provided in the 
results and discussion section. 



where p p , p s , <r p , er s , and p are the local chemical potential 
of polymer and solvent, the deviatoric stress of polymer 
and solvent, and the hydrodynamic pressure, respectively. 
The friction coefficient ( between the two fluids can be 
defined as £ = QtttJs^/^ 3 , where £,(4> p ) is the blob size and 
rj s is the solvent viscosity. The blob size is related to <fi p 
as £ ~ a<p~ m [5], where m = 1 and 3/4 for 0- and good 
solvent, respectively. By minimizing 1Z with respect to v p 
and v Sl the following momentum equations can be derived 



C(f P - v s ) + (/) p Vp + 4> p Vp p - V ■ (t p = 0, (3) 
C(v s - v p ) + (j> s \Jp + <f> s Wp s - V • <r s = 0. (4) 

By combining Eqs.(3) and (4) we obtain 

V-cr-VP-V-7r = 0, (5) 

where the total stress er = er p + er s , P is the modified 
pressure defined by P = p + p s , and rr is the osmotic 
pressure tensor defined as 



V-7T = p (r,<)V^(r,i), 



(0) 



where /i = p p — p s is the difference between the chemical 
potentials. Substituting Eq.(5) into Eq.(4) and eliminat- 
ing P, the following expression for the polymer velocity is 
obtained: 



V • tv - V • er r 



(7) 



2 Theory 



2.1 Two-Fluid Model 

We consider inertialcss fluid motion and free draining poly- 
mers in dilute to semi-dilute polymer solutions. The transient 
evolution of polymer segment volume fraction <f> v is given 
by the equation of continuity: 



In this paper we consider (j) p <^ 1 for dilute and semi-dilute 
solutions, and ^> p V-er s of the above equation is negligible. 
The polymer velocity relative to the total velocity is driven 
by the osmotic pressure and the deviatoric stress of the 
polymer fluid. Substituting Eq.(7) into Eq.(l) yields the 
polymer transport equation [30]: 



dt 



(1) 



where t is time, v p is the velocity of the polymer segments 
in the fluid, 4>p(r, t) = a 3 p p (r, t) with a being the segment 
length and p p the local number density of polymer seg- 
ment. The local volume fraction of solvent is <j) a . Because 
4>p + 0s = 1, the total velocity v(r, t) = 4> p v p + <f) s v s (v s is 
the solvent velocity) satisfies the incompressibility condi- 
tion V • v = 0. The momentum equation of the two-fluid 
model can be derived from the Rayleighian given by Doi 
and Onuki [30]: 



K = 



C(r,t). 



-pV • v 



- ,U p V • {(f> p V p ) - p s V • (</> s V s ) 



i-p : Vv p 



dr, 



(2) 



Ml 
C(0 P ) 



(8) 



This formulation is consistent with Odijk's work [28] ex- 
cept for the (/) p - dependent friction coefficient and the ad- 
ditional V • <T p term. For a solid object with an arbi- 
trary shape, the osmosis-induced force and torque arc / = 
J gn (~ 7T ' n)dS and T = J gn R x (— tv ■ n)dS, respectively, 
where R is a vector to the surface of the object from the 
center of gravity of it. 

Assuming that the polymer stress in the dilute and 



semi-dilute 



< 



overlap fraction 



N 



l-3u 



) polymer 



solution can be expressed as r]'4> p {\/v + (Wv) T ) where i]' is 
a constant, the total velocity follows as v ~ 0(<fip), as seen 



from Eq.(5) and tx ~ O(0p 
Also because p, ~ O(0 P ), so (f> p \7p, p 



and thus V • ov 



0(#). 



O(0p), this implies 



C(V • cr p ) <c O(0 p V/i) and the polymer force V • er p in 
Eq.(8) is indeed negligible [28]. 
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2.2 Chemical Potential of Polymers 

From the self-consistent field theory (SCFT) and the ground 
state approximation (GSA), the Helmholtz free energy of 
the mixture can be expressed as [5] 



F = e 



M 2 + /o( 



>) 



dr, 



(9) 



where e D = k B T/a 3 , C = a 2 /3, f (<f> p ) = (1 - <£ p )ln(l - 
<^p) +X < / ) p(1 — 0p); an< i X is the Flory-Huggins parameter. 
In a good solvent x=0; whereas x=l/2 corresponds to (in- 
solvent conditions. The chemical potential difference (i is 
determined by the functional derivative of F with respect 
to 6 n [281 : 



/j,(r) 



C 



(10) 



Since </> p is much smaller than unity, only the zeroth- and 
first-order terms involving <p p in the chemical potential are 
preserved. The corresponding osmotic stress tensor 7Ty, 
which is scaled by e , can be expressed as 



+c- 



dxi dx-j 



(11) 



where ir n 



fo is the osmotic pressure in 



the bulk. Note that V • ir satisfies Eq.(6). 

Without GSA, the Helmholtz free energy is expressed 



as 



1 



Pl „,A r Q P) 



- W s (f) s 

4> s ln( 



dr 

V0 S 



), (12) 



where A is the number of segment in a single chain, V 
is the volume of the system, and w a and <p a are the di- 
mensionless local interaction field (scaled by k B T) and the 
average volume fraction of polymer segments in the bulk 
for the a-component (a represents polymer or solvent), re- 
spectively. The partition functions Q p and Q s for a single 
chain and a single solvent, respectively, are defined by 



dr / q(r,s)q(r,l — s)ds, 
v Jo 



and 



Q s = 



/ dr cxp(— w s (r)). 
Jv 



(13) 



(14) 



In Eq.(13), q is the statistical weight of the polymer chain 
and can be calculated by 



ds 6 



NV 2 q-Nw p {r)q for 0<s<l (15) 




Fig. 2. A schematic illustration of the convective depletion ef- 
fect on the polymer segment profile (solid line) in a slit between 
two parallel and solvent permeable walls. The dashed line indi- 
cates the equilibrium distribution. The polymer concentration 
vanishes at the walls. 



with the initial condition q(r,0) = 1 and the zero bound- 
ary condition at the solid surface q(r £ dfi) = for an 
arbitrary contour coordinate s. 

In order to evaluate /x p and fi s for a given <fi p (r, t), the 
following iterative procedure can be used [31]: 

i) Find w p (r) and q{r,s) that satisfy Eq.(15) and 

If 1 1 

— / q(r,s)q{r,l-s)ds = —=-<f> p {r,t) (16) 
Q P Jo V(p p 

along with the initial and boundary conditions. 

ii) Evaluate w s (r) that satisfies 

i- exp(-u> s (r)) = -^(r, t). (17) 

Q s V(p s 

iii) Compute the chemical potential difference [i — fi p —fM s 
from the functional derivatives: 

5F 

Mp(f) = j— = -w p (r)+xMr), (18) 
8F 

^ s ( r ) = TT = - w ^ r ) + X4> v {r)- (19) 



iv) Finally, the time evolution of 4> p (r, t) can be calculated 
by Eq.(8) with er p = 0. 



3 Results and Discussion 

We focus on non-equilibrium concentration profiles under 
a uniform flow passing through two parallel and solvent- 
permeable walls (Fig. 2). Both transient and steady state 
results are provided next. In this one-dimensional prob- 
lem, the length is scaled by I = (a 2 /6) 1 / 2 and the time 
scale is defined by r = £ 2 /D, where the diffusion constant 
D = k B T 1 6irr) s ci. The volume fraction of polymer is scaled 
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by the averaged volume fraction <j) p , and the chemical po- 
tential is scaled by k^T /a 3 . Hereafter we use dimension- 
less expressions. The scaled chemical potential difference 
is expressed as 



/' 



<P 



-\7 2 ip + vp 2 , 



(20) 



where v = (1 — 2x)4> P and ip = ((^p) 1 / 2 . Note that the 
ground-state approximation, Eq.(20), for ©-solvent is valid 
only for a narrow slit L/R g < n [5], where R g = a(N/6) 1 I 2 . 
The polymer transport equation (Eq.(8)) thus can be ex- 
pressed as 



at 



Pe 



dip 2 
dx 



B 



-dl-2m]_ 

p dx 



!(2-2m)^M 

dx 



(21) 



where Pe = Ui/D is the Peclet number, and U is the 
uniform velocity. The corresponding boundary conditions 
are <p(0) — and p>(L) = 0. Because the walls are not 
permeable for polymers, the zero flux boundary conditions 
dn/dx\ x —o = and dy/dx\ x= L = are applied. Equation 
(21) leads to the steady state solution: 



H{x) - n(0) 



Pe 



^rn-l)^ dx > (22) 



Under a weak flow condition, Pe < 1, the first-order 
approximation of the steady state solution for ip and y, can 
be written as 

p{x) = ^(j) p {x) ~ ip (x) + Pe <pi(x), (23) 



and 



H(x) ~ Ho{x) + Pe m{x). (24) 



By substituting Eqs.(23) and (24) into Eq.(22), the corre- 
sponding leading and first-order equations become 



d 2 <p 
dx 2 



vtpl + /i o (0)<A> = 0, 



(25) 



and 



^2 +Mo(0) - 3v(p 



<Pi = 



Mi(O) + 0' 



, (26) 



respectively. Both ipo and ipi vanish at the walls. Next we 
consider 0- and good solvent cases under narrow and wide 
slit conditions. 



>< 1 



-e- 




0.4 0.6 

x/R g 

Fig. 3. Steady state polymer segment volume fraction 
(j> p (x,t = oo) of a polymer solution (0 P = 0.1, N = 1000 and 
X = 0.5) in a slit (L/R B = 1, where R g ~ 331) at Pe=0, 0.001, 
0.003, and 0.005. The solid lines are from the numerical solu- 
tions of Eqs.(20) and (21). The symbols (+, x, *, □) indicate 
the analytical approximations, respectively. 



valid and we can use Eq.(20) with v = and Eq.(21) to 
investigate the time evolution of polymer segment concen- 
tration profile. Under steady and weak flow conditions in 
©-solvent (m = 1, v = 0), Eqs.(25) and (26) reduce to 



d 2 <p 
dx 2 



-Mo(0)<A>, 



(27) 



and 



d 2 

^2+^(0) 



-m(Q)<p (x) - <tfp<Po(x) / <p 2 {x')dx'. (28) 



The zeroth-order solution corresponding to the equilib- 
rium state [5] is 



(29) 



with yU o (0) = ir 2 /L 2 , which satisfies the normalization con- 
dition, <f> (x)dx= L. The first order solution reads 

<pi(x) ~ Aain^—x^ 



B< 



L 



1 + 8?^ - -8^1- 



L 



cos(^—x 



3.1 (9-solvent cases 

3.1.1 Narrow slit, L < nR g 

When the distance between two walls L is smaller than 
irRg, i.e., L < irR g , the ground state approximation is 



( 



4 (z*Mz*) ' (30) 



where constant B = (L/tt) 3 p /(16v / 2) and ^(0) in Eq.(28) 
is found to be —<j} p L/2 by the boundary condition <px(L) = 
0. The constant A can be determined by the normalization 
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Fig. 4. (a) Transient evolutions of the segment concentration 
profiles <fi p (x,t), (b) the peak concentration, and (c) the peak 
position. Applied parameters are: P = 0.1, N = 1000, \ ~ 0-5, 
L/R g = 1, and Pe=0.02 at i = t/r > 0+. The profiles shown in 
(a) are at t =0 (no flow), 100, 200, 400 and 1600. At t > 1200, 
the profile approaches a steady state as can be observed in 
panels (b) and (c). 



condition <f>(x)dx = L, and we obtain 



_B 2 
V2 



16 \ 

A = 2ttB + ^= ( 77 - 2tt 2 - — tt 4 J Pe + C(Pe 2 ). (31) 



Accordingly, the peak value and its location can be ex- 
pressed as 

max (Pe) / 14 \ 

=l + W77-2^ + ^7r 4 W + ...(32) 
av (0) V 15 / 



4™ 



and 



l + 2v / 27rSPe + 0(Pe 3 )l. (33) 



To demonstrate how well the first-order approximation 
describes the change of the concentration profile under 
a uniform flow, we consider the steady state case with 
L/Rg = 1 under various values of the Pcclet number and 
compare the approximated profiles with numerical results. 
In Fig. 3, we plot the steady state concentration profiles of 



to 

E a. 

•©- 

aT 
D. 



E cl 

-e- 



CL 



X 

C\J 




0.08 



Fig. 5. (a) Steady state volume fraction profiles of poly- 
mer segments <f> p (x,oo) in the polymer solution (<f> p = 0.1, 
N = 1000, x = 0.5 and L/Rg = 1) at Pe= 0.0, 0.01, 0.02, 
0.04, 0.06, and 0.08. (b, c) Peak polymer segment concentra- 
tion 0™ ax (Pe) and peak position a; ma x(Pe) for various Pe values. 
The filled circles denotes the data obtained by numerical cal- 
culation using Eqs.(20) and (21) and the dotted lines follow 
from the perturbation method results of Eqs.(32) and (33). 



polymer segment under small Peclet number flows, i.e., 
Pe < 0.005 in a polymer solution of N = 1000, <f> p = 0.1 
and x — 0.5. Note that p = 0.1 is slightly below the 
overlap volume fraction 0* ~ 0.11 for polymer chains with 
N = 1000. The analytical results are in good agreement 
with the numerical calculation for Pe < 0.003, but start 
to deviate at Pe = 0.005. Under weak flow conditions 
the first-order solution thus gives accurate results only for 
Pe < 0.003. 

Figure 4 shows the transient behavior of (j) p (x,t) in 
a slit with L/R g = 1 after imposing a fluid flow with 
Pe=0.02 to an equilibrium profile at t = 0. The time evo- 
lutions of the peak value ^>™ ax (£) and its position x max (t) 
are shown in Fig. 4(b) and (c), respectively. The shift of 
the peak position to the steady state position is faster 
than that of the peak concentration, i.e., firstly the peak 
position almost reaches to the steady state position, and 
then the profile of polymer segment concentration be- 
comes sharper. This indicates that a strong convective ef- 
fect applies to the polymer segments and relatively slow 



6 



Takashi Taniguchi et al.: How Flow Changes Polymer Depletion in a Slit 




2 4 6 8 10 




50 100 150 200 250 300 



tx 10 3 /t 

Fig. 6. (a) Transient evolutions of the segment concentration 
profiles 4> p (x,t), (b) the peak concentration, and (c) the peak 
position. Applied parameters are: (j> v = 0.1, N = 1000, \ ~ 0-5, 
L/R g = 10, and Pe=0.002 at t = t/r > 0+. The profiles are 
at t =0 (no flow), 10, 25, 50, 100, 200, and 300 xlO 3 . At 
i = 200 x 10 3 , the profile reaches steady state as seen from (b) 
and (c), and the profiles of t = 200 x 10 s and t = 300 x 10 3 , are 
completely overlapped. Note that the results shown in Fig. 4 
are based on GSA, here the profiles are calculated by using 
the dynamic SCFT scheme. 



relaxation of the polymer distribution across the slit. Next 
we show the steady state concentration profiles p (a;,oo) 
evolve upon increasing the flow rate characterized by the 
range of Pe=0 L 0.01, 0.02, 0.04, 0.06, and 0.08 in a polymer 
solution with <f> p = 0.1, N = 1000 and \ = 0.5 in Fig. 5(a). 
Figure 5(b) and (c) show the peak height </)™ ax (Pe) and 
the peak position x max (Pe) of the concentration profile at 
steady state under various flow strength. The peak height 
^max mcreases quadratically with Pe for small Pe as ex- 
pected in Eq.(32). Above Pe ~ 0.05, </)™ ax increases lin- 
early with Pe. In contrast to the peak height </>™ ax , the 
shift of the peak position for Pe < 0.005 is well described 
by Eq.(33) as seen from Fig.5(b). For Pe > 0.01, the shift 
of peak position to the downstream side is suppressed. 
This is due to the solvent permeable wall and the accu- 
mulation of polymer chains at the downstream side. 




2 4 6 8 10 




=J 0.6 ■ 
1 0.4 ■ 



S 0.2 ■ 

X E 0.0 t ■ ■ ' ' ■ 1 

w 0.0005 0.001 0.0015 0.002 0.0025 0.003 

Pe 

Fig. 7. (a) Steady state volume fraction profiles of polymer 
segments in the polymer solution (<j> p = 0.1, N — 1000, x ~ 0-5 
and L/Rg = 10) at Pe= 0, 1, 3, 5, 10, 20, and 30 xlO" 4 . (b, 
c) Change of the peak concentration <?!>™ ax (Pe) and the peak 
position x max (Pe) for various Pe. The data are obtained by 
using the dynamic SCFT scheme. 

3.1.2 Wide slit, L > nR g 

When the distance L is larger than 7ri? g , we can not rely 
on GSA to evaluate the chemical potential difference n(r). 
Here we utilize the scheme (i)-(iv) based on the self-consistent 
field theory explained in Sec. 2. 2. To demonstrate the scheme, 
we consider the case of L/R g =10. Firstly, we show the 
time evolution of the concentration profiles under Pe =0.002 
in Fig. 6, in which the volume fraction profile of the poly- 
mer solution (x = 0.5) at a quiescent state (Pe=0) co- 
incides with the analytical result of the depletion profile 
near a single wall [32-34], and the peak appears near the 
end of the depletion zone along the downstream side. In 
comparison with the narrow slit case, the shift of the peak 
a^max(i) to the steady state position is much faster than 
that of the peak height (Fig. 6(b) and (c)). Although this 
tendency has been found in the narrow slit case (Fig. 4), 
it is much more enhanced in the wider slit case as seen 
in Fig. 6. Figure 7 shows the steady state concentration 
profiles for Pe=0, 0.1, 0.3, 0.5, 1, 2, and 3xl0~ 3 in the 
polymer solution of <fr p = 0.1, = 1000, \ = 0.5, and slit 
width L/Rg = 10. These are numerical results obtained 
by the dynamic SCFT scheme described in Eq.(8) and 
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Fig. 8. Concentration profiles of s-th segment in the polymer 
solution (<f> p = 0.1, iV = fOOO and \ = 0-5) a t steady states 
under (a) Pe=0 (no flow) and (b) Pe=0.003 in the region 5 < 
x/Rg < 10. The lines are for s = (solid line), 0.05 (dashed), 
0.2 (dash dotted) and 0.5 (dotted). The two insets in (b) are 
the magnified figures of profiles near the downstream-side wall. 



Eqs.(12)-(19). In addition, the peak concentration value 
and its position at steady state are shown as functions of 
Pe. For Pe > lxlO -3 the peak position does not change 
much, but the peak height continues to increase with Pe, 
which indicates that the width of the segment accumu- 
lation becomes narrower while keeping the peak position 
almost the same as Pe increases, showing a strong deple- 
tion effect from the wall. 

The individual segment concentration provides further 
details of the convective effect. From the statistical weight 
q(x,s) obtained by Eq.(15) through the procedures (i)- 
(iv) we can calculate the concentration profile 4>p( s ,x,t) 
for the s-th segment from 



<j) p (s,x,t) 



^q(s,x)q(l-s,x). (34) 



Figure 8 shows the concentration profiles for the s-th seg- 
ment at steady state of s —0.0 (end point), 0.05, 0.20, 
and 0.5 (midpoint) under two conditions: (a) Pe=0 (no 
flow) and (b) Pe=0.003. Only the region 5 < x/R g < 10 
is displayed. The profile for Pe=0 is symmetric about 
x/R g =5, and for Pe=0.003 the concentration is quite small 
for < x/R g < 5. The difference of the distributions 
among polymer segment concentrations is significantly en- 
hanced by the flow. Furthermore, in the case of Pe=0.003, 
the distribution of the end segment (s = 0) is the broad- 
est among others and the distributions for 0.1 < s < 0.5 
is similar to the case of s = 0.5 (midpoint) [2,35,36]. 



3.2 Good-solvent cases 

For good-solvent condition (v > 0) under steady flow, the 
concentration profile can be obtained from Eq.(22) with 




0.4 0.6 0.8 

x/R g 



E Q. 

-©- 
I 



03 
CL 



E Q. 



Q_ 



X 
CM 



6.0 

5.0 

4.0 

3.0 

2.0 

1.0 

0.0 
0.8 

0.6 

0.4 

0.2 

0.0 



.(c) 
















■(d) 


1 1 — 


h- 












0.4 








/>••""" 

1/ 


0.2 
0.0 










0.005 



0.02 0.04 0.06 0.08 0.1 

Pe 



Fig. 9. Results for good solvent conditions, (a) Equilibrium 
(Pe=0) concentration profiles of polymer segments in good 
solvents with \ — —0.5, 0.0, 0.2, 0.3, 0.4, and in ©-solvent 
(X = 0.5) from the bottom to the top lines near the middle 
point, (b, c, d) Steady state behaviors of polymer solutions 
(<f)p = 0.1, N — 1000) in a narrow slit with L/R & = 1 under 
steady flows with various Pe. The polymer segment volume 
fraction profiles (b) in x = are for Pe=0.0, 0.001, 0.01, 0.03, 
0.05, 0.1, 0.2, from left to right. The peak height (c) and the 
peak position (d) for \~ (solid line), 0.3 (dashed line), are 
compared with the ©-solvent (\ = 0.5, m — 1) results, shown 
by the dotted line. As a reference, the case of \ ~ 0-5, but with 
m — 3/4 is shown by the dash-dotted line. 
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rn = 3/4. When Pe=0, we obtain the analytical concentra- 
tion profile as (see Fig. 9(a), details are given in Appendix 
A) 

</> p (x) = [<p m sn(x,k)] , (35) 



where ip m is the maximum value of y/ <f> p (x), sn(cc, k) is the 
Jacobi elliptic integral, and x = xtp m y/v/2k 2 . The con- 
stant k is defined by k 2 = v^/[2/x o (0) — v<p^]. Therefore 
from Eq.(26) we obtain the first-order asymptotic equa- 
tion for the weak flow condition (Pe <C 1): 



— +k 2 + 1 -6fc 2 sn 2 (i,fc) 



dx 2 
L 2 



4K 2 (k) 



sn(i, k) 



Mi(0) 



sn(x' , k)dx 



(36) 



where ipi = <fi/<p m . The constant fii(0) is determined 

by the normalization condition J Q ip 2 (x)dx = L. Unlike 
the (9-solvent case, analytic approximations are extremely 
involved in this case, and thus concentration profiles in 
good solvents are obtained numerically for both weak and 
strong flow conditions. We again consider both narrow- 
and wide-slit cases next. 



3.2.1 Narrow slit, L/R g = 1 

In Fig. 9(b), we plot steady state volume fraction profiles 
(j> p (x) of a polymer solution in a narrow slit with L/R g = 1 
under various Pe. In Fig. 9(c) and (d), both the peak height 
and position are plotted as functions of Pe for various x~ 
values. As a reference the ©-solvent (\ = 0.5, m = 1) 
result given in Fig. 5 is shown by the dotted lines. In ad- 
dition, the case of x = 0.5, but m = 3/4 is shown by 
the dash-dotted line for a comparison. From the compar- 
ison among the cases of x = 0, 0.3, and 0.5 with m=3/4 
in Fig. 9(c), we can see that the peak height of polymers 
concentration in a good solvent become less-sensitive to 
the convective effect with decreasing %. This is because of 
the excluded volume effect in the good solvents. On the 
other hand, the behaviors of the peak height and position 
in the (9-solvent are different from those in good-solvent 
conditions. The different behavior between the good- and 
(9-solvents comes from the p -dependent diffusion coeffi- 
cient in Eq.(21). Namely, the flow in good solvents has a 
larger effective Peclet number by the factor l/^p) 1 / 2 as 
seen from Eq.(21). 



3.2.2 Wide slit, L/R„ 

' I to 



10 



Figure 10 shows the results for the same polymer solution 
in a wide slit with L/R a = 10 and < Pe < 0.01. The 
peak concentration and position are plotted in Figures (b) 
and (c). As seen from (c), the profile is almost flat in the 
quiescent state and by applying a flow the peak position 
jumps from the center of the slit to a downstream position 
near the wall. The polymer segment profile is strongly in- 
fluenced by the applied flow. Because the depletion width 




0.002 0.004 0.006 0.008 0.01 

Pe 

Fig. 10. (a) Steady state profiles of a polymer solution (0 P = 
0.1, N = 1000 and \ = 0) in a wide slit with L/R g = 10 under 
steady flows with Pe=0, 0.5, 1, 3, 5, 7, and 10 xlO" 3 , from 
left to right, (b) the peak height <^>p ax (Pe) and (c) the peak 
position x max (Pe) of 4> p (x,oa) as functions of Pe. The inset 
shows the result in the vicinity of the wall. 



under good solvent conditions is relatively small compared 
with the slit width, the profile near the peak is very sharp. 
This is because the depiction thickness is reduced for good 
solvency by the larger excluded volume effect among the 
polymer chains. 



4 Conclusions 

A theoretical framework is presented based on a combina- 
tion of the two-fluid model and self-consistent field theory 
to describe the polymer segment dynamics as a function of 
a flow field. We investigate a flow-induced concentration 
gradient in one direction only to illustrate the concept and 
demonstrate how the applied flow influences the polymer 
segment concentration profile under Q- and good solvent 
conditions. 

In the case of narrow slits, the ground state approxi- 
mation is valid for Q- and good solvent conditions. Under 
(9-solvent conditions, we obtain analytic expression for the 
steady state volume fraction profile, valid up to the first 
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order in Peclct number Pe. For Pe > 0.003 flow signifi- 
cantly affects the polymer segment concentration profile 
between two walls. For the strong flow cases, we find that 
the peak position with respect to the steady state position 
shifts faster than the peak height of the concentration. 
In contrast to the peak height of the polymer concentra- 
tion, the shift of the peak position in the range of small 
Pe( < 0.005) is well described by the first-order analytic 
approximation. Above Pe ~ 0.005, the shift of peak po- 
sition to the downstream side is suppressed. Under good 
solvent conditions, the concentration profile is now much 
more influenced compared to the ©-solvent case. The de- 
pletion effect in a (9-solvcnt is more significant compared 
to the good solvent case. This is because the depletion 
thickness in a good solvent condition is smaller due to 
the excluded volume effect between the polymer segments 
when the polymer concentration goes beyond the dilute 
regime. 

For wide slit cases, the ground state approximation is 
no longer applicable to the ©-solvent condition. We use 
the numerical dynamic self-consistent field theory in or- 
der to resolve the dynamics for such cases. With time the 
peak position compared to the steady state position shifts 
faster than concentration peak does. The concentration 
profile of segments that differ in ranking number along 
the chains shows that the difference among the polymer 
segment concentration distribution becomes larger when 
applying a flow field. For a flow with Pe=0.003, the dis- 
tribution of the end segment is the broadest as compared 
to all other segments that tend to follow the distribution 
of the center-of-mass segments. For good solvency GSA 
is applicable, and the profile is strongly affected by flow. 
The peak position then gets very close to the wall on the 
downstream side. The model reveals transient relaxation 
and steady state features of the convective depletion dy- 
namics. It will be of interest to validate these findings 
experimentally in the future. 
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A Equilibrium profile in good solvents 

The equilibrium profiles of polymer chains in a good sol- 
vent are described by Eq.(25). Multiplying Eq.(25) by 
dip /dx and then integrate it once, gives 

where C is a constant determined by dip /dx\ x =L/2 = 
and ip (L/2) = <p m at the middle point. Hence C = 
<f% a /j v (0)/2 — Therefore, Eq.(37) becomes 



where (p Q = ip /ip m , x = x(p m ^/v/2k 2 , and k 2 is defined 

by 

(39) 



2 Mo (0)/<- U ' 
Integrating Eq.(38) yields 

<p (x) = sn(x,k), (40) 

where sn(x,k) is the Jacobi elliptic integral. From (p Q = 1 
at x = L/2 we find L = 2K(k) where K (k) is the complete 
elliptic integral defined as 



K(k) = f - = 
Jo J(l-t 2 )(l 



/0 v/(l-t 2 )(l-fc^2) 

and ip m is expressed as 

2K(k) 12k 2 



L V v 



(41) 



(42) 



L /2 

From the normalization condition J <f> a (x)dx = L/2, 
we have 



vL 2 Jo 



sn 2 (i, k)dx = 1. 



(43) 



Based on the properties of elliptic integrals it follows 
K(k)(K(k) - E(k)) = -L 2 , 



(44) 



where E(k) is the second kind complete elliptic integral, 
expressed as 



.•71-/2 



Elk) = f \Jl-k 2 sin 2 Lp dip. 
Jo 

From Eqs.(39) and (42) we obtain 

Mo(0) = — ^(fc 2 + l). 



(45) 



(46) 



(38) 



In summary, for a given v and L wc can determine k 
from Eq.(44), and then evaluate tp m and ^ o (0) from (42) 
and (46), respectively. The equilibrium profiles are zeroth- 
ordcr solutions. 
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